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Abstract 

For the lagrangian L = GlnG where G is the Gauss-Bonnet cur- 
vature scalar we deduce the field equation and solve it in closed form 
for 3-flat Friedman models using a statefinder parametrization. 

Further we show, that among all lagrangians F{G) this L is the 
only one not having the form G r with a real constant r but possessing a 
scale-invariant field equation. This turns out to be one of its analogies 
to /(-R)-theories in 2-dimensional space-time. 

In the appendix, we systematically list several formulas for the 
decomposition of the Riemann tensor in arbitrary dimensions n, which 
are applied in the main deduction for n = 4. 

1 Introduction 



Fourth-order gravity has been a serious alternative to General Relativity since 
1918 already: H. Weyl, see [1], was guided by the idea of the scale invariance 
of the action which requested for an i? 2 -term in its integrand instead of the 
Einstein-Hilbert action integrand R. In fact, the integral / R n ^—gd k x in 
fc-dimensional space-time is scale-invariant just in the case k = 2n, leading 
to n = 2 for the usual space-time dimension k — 4. For details see e.g. the 
reviews [2], [3], [4], and [5], and the books [6] and [7]. For a broader view to 
this topic, and also on the growth of (quantum) perturbations to the today's 
observed large-scale structures by inflation, see the references cited there. 
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Since 1947 it became clearer, that the cosmological evolution can be better 
modeled if both R and i? 2 -terms belong to the action, see C. Gregory [8]. 
In the eighties, the inflationary cosmology has been related to fourth-order 
gravity by Starobinsky [9] , and this paper initiated several follow-up papers, 
e.g. [10] and [11]; generalizations by inclusion of i? 3 -terms and later by a 
general f(R) have been worked out e.g. in [12] and [13]. 

In 1921, R. Bach [14], see [15] for details, initiated a detailed investigation 
of the conformally invariant field equations following from the lagrangian 
CijkiC ljkl . In 1977 it was shown, that a theory with lagrangian of the form 



where units are chosen that light velocity equals 1 and Newton's constant 
equals (IQtt)^ 1 , can be renormalized, see K. Stelle [16]. 

In this context it is often mentioned, that the addition of a multiple of 
the Gauss-Bonnet term, see also [17] and [18], 



to a lagrangian like (1.1) does not alter the field equations, but it leads to a 
surface term which may become essential in the quantization. The relation 
to topology is as follows: the field equations come out by applying contin- 
uous deformations of the metric, but / G^/—gd A x is a topological invariant, 
see [19], [20], [21], and [22]. Lanczos deduced only the 4-dimensional case, 
whereas Lovelock generalized to arbitrary dimensions. His sequence L n starts 
with Li = R, L 2 = G, and each L n leads to a topological invariant in the 
2n-dimensional space or space-time. 

Recently, a lot of papers appeared which contain the Gauss-Bonnet term 
in the action. To circumvent the vanishing of its variational derivative, es- 
sentially three ways have been gone: Models in dimension larger than 4, see 
e.g. [23], [24], [25], [26], [27], and [28], models where G is multiplied by a 
scalar 0, see [29], [30], and models where F(G) instead of G is used in the 
lagrangian with a suitably chosen non-linear function F, see e.g. [31]. Appli- 
cations of theories with Gauss-Bonnet term to cosmology can also be found 



in [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42], [43], and [44]. 



A + R + aR 2 + f3C ijkl C ijkl 



(1.1) 



G = RijkiR^ kl — ARijR^ + _R 2 



(1.2) 



2 



2 Statefinder parametrization 



The metric of a 3-flat Friedman model with synchronized time coordinate t 
reads 

ds 2 = dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) , a(t) > . (2.1) 

We assume that the Taylor development of a(t) exists, the dot denotes d/dt, 
and the Hubble parameter h is defined as usual via 

h(t) = - = a, a = lna. (2.2) 
a 

In what follows, we always exclude a constant function a(t) as it represents 
the trivial Minkowski space-time solution. So, we restrict to functions a(t) 
which have h(t) = at isolated moments of time only; at those moments, 
our exact solutions to be deduced below, have to be matched together. 

A time-inversion leads to a change of the sign of h, so we may assume in 
the following, that always h(t) > 0. Under these circumstances we define for 
any natural number n > 2 

„(n) n n-l m 

z n = , aW = • (2.3) 

d" ' dt n y ' 

This expression is, up to a constant factor, uniquely determined by the con- 
ditions that it is proportional to the n-th time derivative of a with propor- 
tionality factor containing a and a only, is time-reparametrization invariant, 
and is scale-invariant. By obvious reasons it proves useful to define also 
z = z x = 1. 

Another, but equivalent method to define the parameters z n goes as fol- 
lows: it is the only product of with powers of a and a which is dimen- 
sionless in both interpretations of metric (2.1). In the first interpretation of 
(2.1), a is a dimensionless quantity, dimensions are encoded in t, x, y and z. 
In the second interpretation of (2.1), x, y and z are dimensionless quantities, 
dimensions are encoded in t and a; t is being measured in seconds. 

These parameters z n are especially useful if one wants to solve a scale- 
invariant field equation as we are going to do below. As usual, we define 
a field equation to be scale-invariant if for any of its solutions and any 
real constant c, also the homothetically related metric e 2c gij represents a 
solution. 
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Our parameters Z n £1X6 related to the more usual notation, see [45], [46], 
[47], [48], and [49], as follows: 

da . , , 

Z2 = -T7f = —q , z 3 = j , z 4 = -k. (2.4) 
a 1 

Here, q is the deceleration parameter, j the jerk parameter, and k the kerk 
parameter. The notion statefinder parameter refers to the pair (j, s), where 
s is defined for q ^ 1/2 via 

1 ' (2.5) 



3(2-1/2) 



Next, we give some relations between the parameters z n and the Hubble 
parameter h: Solving eq. (2.3) for we get 



flW = %^ . (2.6) 



a 



n-l 



The temporal derivative of eq. (2.6) has the same l.h.s. as eq. (2.6) with n 
replaced by n + 1. Equating the related r.h. sides we get 

z n +i = z n /h + z n (nz 2 + 1 - n) . (2.7) 

From eqs. (2.2) and (2.4) we easily deduce z 2 in dependence of h: 

z 2 = -q = l-^( l -]=l + h/h\ (2.8) 



dt \h 

and for n > 3 we can iteratively deduce z n in dependence of h with eq. (2.7), 
the next two terms being 

z 3 = j = ^ + z 2 {2z 2 - 1) = 1 + 3h/h 2 + h/h 3 (2.9) 

and 

6h 4/i 3^ 1 d 3 h 

So we get the relation to the other set of dimensionless constants, see eq. 
(3.12) of [2] 

£ -?Ll. h -P-\ (2.10) 

p dtp y ' 

This leads to Z\ = e = 1, z 2 = 1 + £i, and ^ 3 = 1 + 3£i + e 2 . 
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If we take the logarithmic cosmic scale factor a (see eq. (2.2)) as new 
time coordinate we can rewrite eq. (2.7) as 

dz n 

z n +i = -r- + z n (nz 2 + 1 - n) . (2.11) 
da 

After some reformulation we also get as metric 

ds l = _ ( dx l + dy 2 + d ^ (212) 

with the scale invariant parameters being 

1 dh / \ 

'=-'-rp (2 - 13) 

from eq. (2.8), and from eqs. (2.9.) and (2.10) 

J = 2g 2 + g-^ (2.14) 
da 

and 

k = 3jq + 2j-^. (2.15) 
We will apply these formulas in section 4. 



3 Gauss-Bonnet lagrangian 

For 2-dimensional space-times, lagrangians of the type f(R) have been dis- 
cussed e.g. in [51], [52], [53], and [54]. In [51], the lagrangian 

f{R) = R k+1 (3.1) 

was shown to lead to non-trivial classical results even in the limit k — > 0. In 
[53] , this limit was shown to produce the same field equation as the lagrangian 

f(R)=R-\nR. (3.2) 

This property is related to the fact that / R^fgd 2 x is a topological invariant 
related to the genus of the space. 

Similarly, in [54], the integrand R was kept constant, but instead, the 
dimension of space-time was formally defined as 2 + e, and the limit e — > 
was discussed. 
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We now want to transfer this idea to the set of 4- dimensional space-times. 
To this end we consider a general function F(G) with G from eq. (1.2) as 
integrand of the action. The full field equations are given e.g. in eq. (5) of 
[35] using the notation F G = dF{G)/dG, 

= \g ij F{G) - 2F G RR ij + AF G R l k R kj - 2F G R lklm R j klm 

-AF G R lkl3 R H + 2RFg - 2g ij RDF G - AR ik F^. k 
-AR jk F£. k + m ij UF G + Ag^R kl F G . kl - AR^ l F G , kl . (3.3) 

The integral I G — J Gyf^gdfix is a topological invariant related to the Euler 
characteristic. Therefore, the function F(G) = G leads to the field equa- 
tion reading = trivially fulfilled by all metrics g^, i.e., every space-time 
represents a stationary point of the action I G . 
Next, we consider the action 

1 = J F(G)v^rf 4 x (3.4) 

and ask for its properties if a scale transformation is applied to the metric. 
More exactly: What happens with I eq. (3.4) if we replace g^ by its homo- 
thetically equivalent metric e 2 " 1 g^ where 7 is an arbitrary constant? If / does 
not change at all by this transformation, then we call / scale-invariant. G goes 
over to e~ 4l G by this transformation, and so, obviously, only F(G) = c 2 ■ G 
with a constant c 2 leads to a scale-invariant action J. 

A less trivial question is the following one: Under which conditions, the 
action / eq. (3.4) is almost scale- invariant, i.e., scale-invariant up to adding 
a multiple of I G 7 In other words: Which functions F(G) have the property 
that replacing g^ by e 2n ' g^ in eq. (3.4) leads to the action I + kj ■ I G with 
constant A; 7 ? The answer is: Besides the case already discussed above, 

F(G) = Cl -G-\nG + c 2 -G (3.5) 

with constants c\ 7^ and c 2 is the complete set of solutions. 1 A word 
to dimensions: the argument of the logarithm should be dimensionless, so, 
instead of InG we should have written ln(G/G ). However, a change of the 
value G can be compensated by a redefinition of the constant c 2 . As the 
term with c 2 does not contribute to the field equation, we may put it to 

1 For negative values of G, the term lnG should be replaced by In \G\. The singularity 
at G — > is a mild one and in the models of our interest, \G\ is positive anyhow. 
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zero classically. Dividing everything by c\ we finally get the only interesting 
remaining almost scale-invariant case to be 

F[G) = G ■ In G . (3.6) 

If we insert eq. (3.6) into eq. (3.3), the following simpler field equation 
appears: 

= -g^G ■ In G - 2(RR ij - 2R l k R kj + R tklm R\ lm 

+2R mj R kl ) ■ (1 + InG) + 2R{\aG)' ij - 2g ij Rn(\nG) - AR ik (\nG)% 

-AW k {\nG)% + m ij U(hiG) + 4g i3 R kl {lnG). M - AR ikjl (\nG). kl . (3.7) 

Due to its importance it seems justified to deduce this case by another way: 
Take a small positive parameter e and define 

F e (G) = - ■ (G 1+t - G) (3.8) 

which leads to the same vacuum equation as the lagrangian G 1+e . Then the 
limit e — > in eq. (3.8) exactly leads to (3.6). Sketch of the proof: Put 
G = e x , then 

G e = e tx ^l + ex = l + e\nG. 



4 Exact Friedman models 

We apply the notation of section 2, especially metric (2.1) with Hubble pa- 
rameter (2.2) etc. If we start with a(t) = t n with positive values n and 
t, we get h = n/t, a = nlnt, q = (1 — n)/n, j = (n — l)(n — 2)/n 2 , 
k = —(n - l)(n - 2)(n - 3)/n 3 , t = e a/n , and h(a) = n ■ e~ a/n . The metric 
can then also be written as 

ds 2 = ^ e ^/n _ e 2a ( dx 2 + dy 2 + ^ 

This leads to the de Sitter space-time as n — > oo where q — — 1, j — 1, 
fc = —1, and s = 0. Within Einstein's theory and with pressureless matter 
of density p, the deceleration q is related to the critical density p c necessary 
to close the universe via 2q = p/p c - 

Using eqs. (6) and (7) of [35], or using eq. (3) of [38] we get 

R = 6 (h + 2h 2 ) , G = 24 (hh 2 + h A ) (4.2) 
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and the vacuum equation following from the action (3.4) as 

= G ■ F G - F(G) - 24G • h 3 ■ F GG , (4.3) 

where F G = dF/dG and F GG = dF G /dG. In comparison with the full field 
equation (3.3), this is a surprisingly simple equation. We test the previously 
discussed property as follows: adding C2 • G to this F, the set of solutions to 
eq. (4.3) will not change. For non- vanishing F GG , i.e., a non-linear function 
F(G), eq. (4.3) is of third order in the metric, as it represents the constraint 
equation to the full fourth-order field equation. 

Now we insert the example F(G) = GlnG of eq. (3.6) into eq. (4.3) 
and get via F G = 1 + lnG and F GG = 1/G and after multiplication with 
G = -24/i 4 • q 

= G 2 - 2AG ■ h 3 . (4.4) 

The singular case G = needs an extra consideration: Looking at eq. (4.2) 
this leads to h — —h 2 , as the case h — was already excluded earlier. This 
behaviour can be written in the original form (2.1) by a(t) = t, i.e. the 
deceleration vanishes identically: q = 0. 

Now we look for the remaining solutions of eq. (4.4), i.e., those with 
G 7^ 0. To this end we insert eqs. (2.8) and (4.2) into eq. (4.4). The result 
is the second order equation for h 

(hh 2 + h 4 f = h 3 ■ ^(hh 2 + h A ) (4.5) 

reducing via f = see eq. (2.2), to the following first-order equation for 
the deceleration parameter q 

^- = 4q + 3q 2 . (4.6) 
da 

The fact, that eq. (4.6) does not contain the Hubble parameter is a 
consequence of the scale-invariance of the field equation. By the way, eqs. 
(2.14) and (4.6) can be combined to q 2 + 3q + j = characterizing this field 
equation. 

The other solution with constant value q is q = —4/3. Using eq. (2.8) we 
get 3h = h 2 , i.e. h = — 3/t and finally a(t) = 1/t 3 . These two solutions with 
constant q, i.e. a(t) = t and a(t) = 1/t 3 , represent themselves self-similar 
space-times: Multiplying the metric of space-time with a constant factor can 
be compensated by a time-translation. 



8 



Let us finally come to the case on non-constant g in eq. (4.6). Considering 
solutions as same, if they are related by a scale-transformation, exactly three 
solutions remain, characterized by 



«<«> = -3 + sTF=' -t <<!< ° (47) 



and 



where a may take all real values in eq. (4.7), but eq. (4.8) is not defined for 
a = and represents one solution for a > 0, i.e. q < —4/3 and another one 
for a < 0, i.e. q > 0. Coming back to a relation for the scale factor, we get 

/ \ da 4 

= =-ii±w- (49) 

One of the two remaining quadratures can still be done in explicit form via 

:(ln h 
da 

C 

a 



eq. (2.13), i.e. = -1 - q, leading to 



%) = -V±ir V3 (4-10) 



with a positive constant c. The final step to get the function a(t) is then via 
the integral 

r°(*) , 

\a 4 ±l\- 1/3 da = c-t. (4.11) 



/o(0) 

5 Conclusion 

For the lagrangian L = GlnG where G, see eq. (1.2), is the Gauss-Bonnet 
curvature scalar we deduced the field equation and solved it completely up 
to one final quadrature eq. (4.11) in closed form for 3-flat Friedman models 
using a statefinder parametrization. Further we have shown, that among all 
lagrangians F(G) this L is the only one not having the form G r with a real 
constant r but possessing a scale-invariant field equation. This turns out to 
be one of its analogies to /(i?) -theories in 2-dimensional space-time. 

Recently, several other modifications of Einstein gravity have been dis- 
cussed, see e.g. [55] for a non-local one, here we propose with the arguments 
given above, as gravitational lagrangian 

L g = A + R + aR 2 + pC ijkl C ijkl + 7 G In G (5.1) 

being worth considered in more details than done up to now. 
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6 Appendix 



Here we present some decompositions of the Riemann tensor from the geo- 
metric point of view which are implicitly used in the text above, and which 
may have some interest in themselves and have other applications, too. 

In four dimensions, the Riemann tensor Rijki possesses 4 4 = 256 real 
components. By use of the known symmetries, this figure reduces to 20, but 
this twenty-dimensional space is even harder to imagine. Example: To work 
with the field equation (3.3) it is necessary to know, that in four dimensions, 

(~iiklrri(~< ri /~igklm(~< 

^ ^jklm — ^ Oj O ^gklm 

and how this can be used for evaluating analogous terms with the Riemann 
tensor. 

Below, we will present four different possibilities how to arrange this set 
of components to get a better understandable system. 

The Riemann tensor R^m of a space-time of dimension n > 3 can be 
decomposed according to several different criteria: 

1. The usual one into the Weyl tensor Cijki plus a term containing the Ricci 
tensor R^ plus a term containing the Riemann curvature scalar R. 

2. Two trace- less parts plus the trace. 

3. The Weyl tensor plus only one additional term. 

4. Two divergence-free parts plus the trace. 

We use the following two properties of the Riemann tensor 

Rijki = —Rijik , Rijki = Rkiij ■ (6.1) 

The Ricci tensor is the trace of the Riemann tensor: R^ = g M R ik ji, where 
gki denotes the metric of the space-time, and the Riemann curvature scalar 
is the trace of the Ricci tensor R = g kl R k i- The sign conventions are defined 
such that in Euclidean signature, the curvature scalar of the standard sphere 
is positive. 

For any symmetric tensor Hij we define another tensor H*j kl via 

H*jki — Hi k gji + Hjigi k — Hugj k — H jk gu . (6.2) 

Then the tensor H*j kl automatically fulfils the identities eq. (6.1). For the 
special case = we get the simplified form 

9ijki = 2 9ik9ji - 2gu9jk ■ (6-3) 
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6.1 The usual decomposition 

The Weyl tensor C^ki is the trace-less part of the Riemann tensor, i.e. 
g lk Cijki = 0. It vanishes identically for n — 3. Using the notation of eqs. 
(6.2) and (6.3) we make the ansatz 

Rijki = C ijkl + aR* jkl + BRg* ]kl . (6.4) 

Then the coefficients a and B have to be specified such that the trace-lessness 
condition for the Weyl tensor becomes an identity. This condition determines 
the coefficients a and 8 uniquely, and the result is: 

a = and 8 = — — — . (6.5) 

n-2 H 2(n-l)(ra-2) v ; 

Thus, we get the usual formula 

Rijkl = Cij k i H (RikQjl + RjlQik — RiWjk — RjkQil) 

1 

R (dikdjl — 9u9jk) ■ 



(n-l)(n-2)' 

6.2 The decomposition using trace-less parts 

In distinction to the previous subsection, we now perform a more consequent 
decomposition into trace and trace-less parts. To this end we define as 
the trace-less part of the Ricci tensor, i.e. g^Sij = with Sij = Rij + nRg^ 
possessing the unique solution k = — 1/n, i.e., = R^ — Rgij/n. Then the 
analogous equation to eq. (6.4) is 

Rijki = Cij k i + lS* jM + 5R g* jkl . (6.6) 

This becomes a correct identity if and only if 

7 = and 5 = — : . (6.7) 

' n-2 2n{n - 1) v ; 



So we get 



Rijki — Cijki H (S ik gji + Sjig ik — Sugj k — Sj k gu) 

+ —, — ~jT^ (dik9jl - 9il9jk) ■ 



n 



6.3 Decomposition into two parts 

Let us define a tensor Ly = R^ + (R such that a parameter e exists which 
makes 

Rijki = Cij k i + eL* jkl (6.8) 
becoming a true identity. It turns out that this is possible if and only if 

C = . ~ l r and e = — !— - . (6.9) 
s 2(n-l) n-2 V ; 

Thus, we can write = R^ — 2 (n-i) R 9ij an d 

Rijki = Cij k i H — — (Li k gji + Ljig ik — Lugj k — L jk gn) . (6.10) 

6.4 Decomposition into divergence-free parts 

Now, besides the identities from eq. (6.1), we also use identities involving the 
covariant derivatives, denoted by a semicolon, of the Riemann tensor. The 
Bianchi identity reads 

Rijkl;m Rijlm;k Rijmk;l . 

Its trace can be obtained by transvection with g lk and reads 

Rjl-m + R l jlm;i ~ R jm;l = . (6.11) 

It should be mentioned, that the transvection with respect to other pairs of 
indices does not lead to further identities. The Einstein tensor is defined 
as Eij = R^ + XR gij , where A has to be chosen such that the Einstein tensor 
is divergence- free, i.e., E t j. i = 0. Using the trace of eq. (6.11) (again, there 
is essentially only one such trace), namely 2R t l . i — R.j = 0, we uniquely get 
A = —1/2, i.e. the Einstein tensor is E^ = R^ — ^Rg^. 
With the ansatz 

Rijki = W ijkl + r)E* jkl + 9R g*j k i (6.12) 

it holds: The coefficients 77 and are uniquely determined by the requirements 
that eq. (6.12) is an identity, and the divergence of the tensor W ijk i vanishes: 
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W l j k i ;i = 0. We get uniquely the following values of the constants: rj — 1 and 
9 = \. Then 

Rijki = Wijki + E ik gji + Ejig ik — E a gj k — Ej k g a 

+ T^R {9ik9ji - 9u9jk) (6. 13) 

defines a decomposition of the Riemann curvature tensor into the divergence- 
free tensors Wij k i, E^, and the scalar R. 

It should be mentioned, that for every n > 2, the four tensors Rij, Sij, 
Lij, and E^ represent four different tensors. And it is a remarkable fact, that 
the coefficients in Eij and in eq. (6.13) do not depend on the dimension n. 
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